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Abstract 
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I. INTRODUCTION 



While the late-time acceleration of the expansion of our universe seems to be a robust 
conclusion reached by several independent cosmological observations [H, Q, 0, 0, [sl , its origin 
remains an open question [gI, 0]. Various explanations for late-time acceleration have been 
proposed ranging from modifications of the laws of physics to questioning the Copernican 
principle jl, U 10|. If this principle holds we are lead to the conclusion that our cosmological 
model requires the introduction of new physical degrees of freedom, often referred to as the 
"dark sector". In this case, it is necessary to specify the interactions of these new fields, both 
their self interactions as well as their interactions with matter. There are several candidates 
for this new matter, the simplest being quintessence models llj, or a modification of general 
relativity (see Ref. [7(] for a discussion on this distinction). Let us emphasize that the simplest 
model of all requires only the introduction of a single constant, the cosmological constant, 
and is at the moment in agreement with all available data j^, l3 • 

Among all the theories of gravitation beyond general relativity, scalar-tensor theo- 



ries [13|, [ij] are a well motivated extension and may also be related to quintessence since they 
include one or more scalar fields universally coupled to matter. These theories are described 
by two functions, a potential describing the scalar self interactions, and a coupling function 
describing interaction of the scalar field to matter. Indeed, these theories depart from gen- 
eral relativity and are subject to many experimental constraints, particularly from the Solar 
system [l5|. These deviations can however be larger on cosmological scales because, as l ong 



as the coupling function has a minimum, such theories can be dynamically attracted Ig 
toward general relativity during the evolution of the universe. More generally, when consid- 
ering theories with a quintessence-like potential, it was shown [13, 3] that the theory can 
be attracted toward general relativity while the scalar field is driven toward the quintessence 
tracking solution, hence leading to models of extended quintessence that were advocated to 
be very generic for quintessence models [l^, 20|. In that context, the deviations from gen- 

the cosmic 



eral relativity can be constrained from big bang nucleosynthesis (BBN) |2ll . 1221 
microwave background (CMB) [2^, weak lensing and the production of gravitational 
waves j25|. Interestingly, they can also lead 26] to an equation of state with p/p < —1 for 
the dark energy. 

It is indeed also possible that the new scalar field(s) do not couple universally to matter. 



as considered in Refs. |19l . |20[ |. In this case, one expects a violation of the equivalence 
principle and a variation of the fundamental constants which can used to constrain the 
equation of state of the dark energy 
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constrained 
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However, such couplings are in general severely 
with the exception of some particular cases such as the coupling to 
neutrinos [30| or dark matter. 

When the scalar fields couple differently to the dark sector, there is almost no constraint 
from the equivalence principle on the interactions between dark matter (DM) and standard 
(or visible) matter (see however Ref. j3l| for a possible test of the universality of free fall 
between DM and normal matter). Such a possibility was initially investigated in Ref. 

33( 1 . These 
and have been 



in the particular case of Brans-Dicke theories and was recently revisited in Refs 
models can, in principle explain, the coincidence problem as well [3. 



35, 36 



argued to naturally appear when the quintessence model is embedded in a supersymmetric 
construction (37j . 

In this article, we focus on big bang nucleosynthesis and build on our previous analy- 



sis |22l . |38| that focused on the constraints that can be imposed on scalar-tensor theories of 
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gravity with a massless dilaton. We extend the results of Ref. 22| in two directions: (1) by 
considering extended quintessence models, i.e. a general scalar-tensor theory with a poten- 
tial, and (2) by extending our formalism to models with a non-universal coupling between 
the dark matter sector and the visible sector. The theories considered in this article are first 
detailed in § [IT] and both cases are studied respectively in §|in]and §[lYl Most technicalities, 
such as the BBN data and Solar system constraints, are gathered in appendices. 



II. SCALAR-TENSOR THEORIES OF GRAVITY 
A. Definition of the model 

In this article, we consider models in which gravity is described by a scalar-tensor theory 
and in which the strength of the scalar coupling may differ between CDM and ordinary 
matter. It follows that in the Einstein frame the action takes the form 

/i4 
TTT^v^ [R* - 2grd,^,d,^, - AViv,)] 

+Sv[AUv*)9;u:i^v] + SD[AUv*)9;u;i^D]. (2.1) 

The action contains three arbitrary functions: the potential V, the couplings to ordinary 
matter. Ay {V for visible sector), and the coupling to CDM, Ad {D for dark sector). We 
will also consider cases where both coupling functions are set equal to each other. is the 
bare gravitational constant from which we define = SttG*. 

The Jordan frame is typically defined as the frame for which the metric is minimally 
coupled to matter. When Ay = Ad, there is a unique definition of the Jordan frame. With 
two coupling functions, we can define two different Jordan frames. Since standard clocks 
and rods are made of standard matter, we define the matter Jordan frame (MJF) to be 
that frame in which the standard matter fields are minimally coupled to the metric. The 
MJF metric will thus define length and time as measured by laboratory apparatus so that 
all observations (time, redshift,. . . ) have their standard interpretation in this frame. This 
is also the frame in which the nuclear reaction rates take their standard form, since they do 
not involve CDM. In this frame the stress-energy tensor of visible matter is conserved. The 
action then takes the form 

+Sv[9f,u; i^v] + SD[B'^{^*)9tMu; tpo] (2.2) 
after performing the conformal transformation 

g;^ = n^)9,u (2.3) 



with 



dv^^V 3 



dlnF( (^) 

dip 



n 2 



+ 2F(^) ^'-"^ 



dip J 4 

Av{v.) = F~'/\ip) (2.5) 

2V{^.) = U{v)F-\^) (2.6) 

B{ip,) = Ad/Av. (2.7) 
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Wc will denote all Einstein-frame quantities with a star (*). The Einstein frame has the 
advantage of diagonalized the kinetic terms for the spin-0 and spin-2 degrees of freedom so 
that the degrees of freedom of the theory are more easily discussed. The strength of the 
couphngs of the scalar field to the matter and CDM fields is characterized by 

, , din , „, 

ai{ip,) = — (2.8) 

d(^* 

and we also define 

d (y • 

(3,M ^ (2.9) 
d</7* 

with i — V,D. 



B. Cosmological background equations 



1. Space-time metrics 



We consider a Friedmann-Leniaitre universe with metric in the MJF 

ds^ = -df + R^{t)^ijdx'dx^ (2.10) 

where ■jij is the spatial metric and R the scale factor. In the Einstein frame, the metric 

takes the same form but with a scale factor and a time coordinate, which are related 
to their MJF counterparts by 

R = Av{ip*)R*, dt ^ Av{ip*)dt^ (2.11) 

and the redshifts are related by 

l + z= + (2.12) 



2. Equations of motion in the Einstein frame 



The Priedmann equations in this frame take a form similar to those of general relativity 
with a minimally coupled scalar field and a fluid, 

3 (^H^, + ^) = SttG^p, + iPl + 2Viip,) (2.13) 
3 d^R 

* = 47rG,{p, + 3P,) + 2ijl - 2V{(p,) (2.14) 



where we have introduced = dlnR^/dt^ and 

-0* = d</7*/di*. (2-15) 

is a positive, null or negative constant specifying the curvature of the spatial sections, 
and P* are the energy density and pressure of the total fluid. 
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In this model, the total energy can be split between the dark sector {pn) and the visible 
sector (py), 

P = Pv + Pd- 

We assume that the dark sector component is described by a constant equation of state wd, 
which is set to wo = for CDM. The visible sector can be further decomposed between a 
pressureless matter component Pmat = and radiation (Prad = -^Prad), so that 

PV = Pmat + Prad- 

We stress that Orad = «mat = by construction, in order to satisfy the universality of free 
fall. In the Einstein frame, the evolution of these energy densities are dictated by the matter 
conservation equation, 

+ 3if*(P« + Pi*) = "i(v^*)(p« - 3Pi*)?/'*. (2-16) 

The r.h.s. term proportional to ip^: is due to the scalar interaction. In our case it reduces to 
0, aypmat*"^* and aDPD*ip* for the radiation, standard matter and CDM. 
The evolution of the scalar field is a Klein-Gordon equation 

+ ?>H^il)^ = --^ - 47rG'* y^ai(v9*)(pi* - 3Pj*) , (2-17) 
dr* acp^, ^—^ 

where the sum is taken over i = rad, mat, D. The coupling term reduces to 

-47rG'* ^«i(v9*)(pi* - 3Pi*) = -AnG^ [aoPD* + «yPmat*] • 



3. Implementing BBN 

The nuclear reaction network takes its standard form in the matter Jordan frame, where 
the Lagrangian for the visible sector is not affected by the existence of the scalar field so 
that the cross-sections etc. . . also take their standard form. 

To compute the light elements abundances during BBN, one only needs to know the 
expansion rate history, H{z), from deep in the radiation era up to today. It is thus convenient 
to express the Hubble parameter in the MJF in terms of the one in the Einstein frame, using 
Eq. dm]), as 

AyH = [H, + avi^M ■ (2-18) 
Eq. (12.181) can also be expressed in the simple form 

(2.19) 

where p is the number of e-folds of expansion in Einstein frame, 

p=-\n{l + z,). (2.20) 

Our numerical strategy follows exactly the one developed in Ref. [13]. We solve the Ein- 
stein equation in the Einstein frame (see appendix|X]) and then deduce H{z) from Eq. (I2.18P 
to compute the production of the light elements. 



AyH = 



1 + av{^* 



|d^ 
dp 
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4- Evolution the energy densities 



The system of equations flAlHATp are expressed in terms of the energy densities of the 
various components. Interestingly, these can be derived analytically. In the MJF, the 
densities are related to their Einstein frame counterpart by 

= A'yp, = A\rP. (2.21) 

For matter in the visible sector, Eq. (12.161) takes the simple form 

pv + 3i/(py + Py) = (2.22) 

and is trivially integrated to get 

py. = py.o(l + ^)'^'+"''^ , (2.23) 

for any component of the visible sector (i.e. i = rad, mat) with a constant equation of state 
and where z is the redshift defined by 1 + z = Rq / R. Using Eqs. (12.11112.121) . we obtain the 
evolution of the energy density in the Einstein frame 

/ A \ 4-3(i+t«,) 

Pv.. = Pvio* -r^ (1 + z,f^'^^'\ (2.24) 



^Vio = — WTTT^Pvio- (2.25) 

O-Hn 



where pio* = ^voPvio- We can also define the associated density parameter by 

In the dark sector, the conservation equation has a source term in the MJF, 

Pd + 3H{pD + Pd) = Ay\aD - av){pD - ^Pd)4^*- (2.26) 

This equation can be obtained by plugging Eq. fl2.2ip into Eq. (12.161) foi i = D and using 
the relation (12.111) to shift to derivative with respect to t*. Interestingly, we can find the 
solution of this equation analytically. In the CDM Jordan frame, defined by t and scale 
factor R, CDM is not coupled to the scalar field so that pz) oc p-^ci+wo)^ QQj^g back to the 
Einstein frame, this implies that 

PD. OC i?;3(l+.,)^4-3(l+.,) 

and then in MJF, 

PD = Pm(4^^\ (l + .)^(^+-). (2.27) 



Indeed, if the coupling of the scalar field is universal, the dark matter and normal matter 
satisfy the same evolution equation. In a general model where Ay ^ Ao-, the two matter 
components evolve differently. This is related to the fact that both forms of matter do not 
experience the same coupling to the scalar field. 

We deduce that the ratio between the dark matter and baryonic matter components is a 
priori not constant during the evolution of the universe. It is given by 

pD* pD PDO AvqAd, , „ P(^*) , . 

- 1 7-[f*) = , [Z.ZH) 

Pm* Pm PmO ^DO -DQ 
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with So = Pool Pmo and we define So = Eq/Bq. This is a main feature of this class of models, 
which distinguishes it from universal scalar-tensor theories. 

Finally, the energy density of the scalar field in the Einstein frame is given from Eq. fl2.13p 

by 



dip. 



so that its energy density in the matter Jordan frame is 



(2.29) 



(2.30) 



III. AN EXTENDED QUINTESSENCE MODEL 

Before we move on to the most general type of models with two coupling functions, 
we consider first an extension of the models considered in Ref. 2^ with self-interaction 



potentials, but with Ay = A^. Our first examples are thus models of extended quintessence 
in which the scalar field accounts for the late time acceleration of the Universe while also 
being responsible for a scalar interaction. Given the framework we have set up in the previous 
section, we can now examine the consequences on BBN of several choices of quintessence 
models. 

As a first example, we consider an extended quintessence model of the runaway type with 
a potential and coupling function given by 

V^ = MVr, \nA = Cif:\ (3.1) 

where a, b and C are all positive and the mass scale M <^ hT^I"^ must be tuned so that 
~ today. It follows that 

a = ~Cb^:'-\ (3.2) 

When C = 0, this model is a standard quintessence model with an inverse power law 
potential. The dynamics of this model has been investigated in depth. In particular, there 
exist scaling solutions with attractors of the dynamics of the scalar field which behaves as a 
perfect fluid with an effective equation of state 

wa~2 

= (3.3) 

where w is the equation of state of the fluid that dominates the matter content of the 
universe. Figure [1] depicts the evolution of the energy densities as a function of the redshift 
in this case and shows the quintessence attractor mechanism typical for inverse power law 
quintessence potentials when gravity is described by general relativity (i.e. C = 0). 

When C 7^ 0, and as discussed in detail in Ref. [2^, we expect two main effects. First 
during BBN, the field dynamics is modified during electron-positron annihilation and then 
the late time dynamics is modified because of the coupling to non-relativistic matter. Let 
us investigate them by assuming that = 0. 

As detailed in § III. A. 2 of Ref. |i22i], the dynamics during electron-positron annihilation 
is governed by the equation 

■5 2 

V* + nV'* = -a(v5*)Se, 
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FIG. 1: Evolution of the energy densities of the matter, radiation and scalar field as a function of 
the redshift for a model with potential (j3.ip with a = 6, C = and M = 4 x 10~^^ for different 
initial conditions, (/>in = 2 X 10-l^ 6 X 10-1^, and 2 10"^^ The field initially enters a kinetic regime 
so that oc (1 + z)^ and then reaches a slow-roll regime until its energy density is of the same 
order of magnitude as that of the matter and radiation. At this stage, it enters a tracking regime 



where its effective equation of state is w^, = and 



-1/4 during the radiation and matter era 



respectively. It follows that the scalar field always ends up dominating the matter content of the 
universe at late times and gives rise to an accelerated phase of expansion. 



where Se = (pe — SPg) / prad and prime denotes a derivative with respect to p. In a non-kinetic 
regime, this equation can be approximated by 

3 



As in the case of a quadratic coupling 13] , the field is frozen to some initial value y^in* during 
the radiation era |39[] prior to the electron-positron annihilation. Then the source term Eg 
acts temporarily inducing the evolution of ip^, which then settles to another constant value 
yjout* in the radiation era. Because In A does not have a minimum, the field will shift towards 
larger values so that v^out* > V^in*- 

During the matter era, the evolution equation reduces to 

Since the minimum of the coupling function is at infinity, it follows that the scalar field is 
simply attracted toward this value and does not undergo damped oscillations as in the case 
of a quadratic coupling. In particular, in the slow-roll regime, it behaves as 



Cb 

l + (2 + &)^(p-Peq) 



1 

2+6 
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where peq is the value of p at the matter-radiation equahty. 

In the general case where both the potential and the coupling are effective, the solution 
is first attracted toward the quintessence scaling solution during the radiation era during 
which the coupling is not efficient. The various mass thresholds further drive ip^ toward 
larger values so that it reaches the quintessence tracking solution more rapidly. Then, 
during the matter era, the evolution is driven by the potential and coupling which both 
drive (f^ to infinity so that the scalar-tensor theory is attracted toward general relativity. 
Figure [2] shows the magnitude of these effects. 




FIG. 2: (Left): Evolution of the energy densities of the matter (blue line), radiation (red line) and 
scalar field as a function of the redshift for a model with potential (j3.ip with cl){z = 10^^) = 1, 
a = b = 6 and C = (solid black line) or C = 0.2 (dash dotted green line). (Right): Evolution of 
6 for the same conditions. 



Quintessence fields coupled to the Maxwell term, F"^, in the action, will induce a variation 
of the fine structure constant. There are a few models kO, 41 



which can simultaneously 

account for the purported variation in quasar absorption systems |42^ and the Oklo and 



meteoritic constraints [43|, |44, |45| . Here we consider the model of Ref . [40] with the potential 



(3.4) 



which is representative of a class of models which possess a minimum in the self-interaction 



potential and is closely related to the one studied in Refs. [H, 13]. The potential is nor- 
malized by taking Vq to correspond to the present vacuum energy density. The coupling 
function in this model is also specified and related to the potential 



b+V{^,)/Vo 
l + b 



(3.5) 



where 1 + 6 > 0. Its main feature is to have a common minimum with that of the potential, 
and allows the system to recover general relativity when the scalar field is at the minimum. 
This is in the spirit of the least coupling principle jTsj whose goal is to suppress the scalar 
interaction when the field has reached its minimum. Note that in Ref. ^3] where the scalar 
field was also responsible for a variation of the fine structure constant, constraints were 
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derived such that n needed to be very smalL In the current context, the field is universally 
coupled so that it is responsible only for a variation of the gravitational constant, which is 
much less constrained than the time variation of the fine structure constant jlsl, 2^ . 



For small values of the field ip^, A can be approximated by 1 + nXipl/ 2(1 + b), so that we 
can write 

a ^ /3 ^ (3.6) 

1 + 61 + 6 ^ ' 

For /5 = 0, we show the evolution of y?* as a function of redshift in Figure [3^. We see that 
for most of the evolution yj* is constant largely due to the very small value of Vq. At late 
times, the field evolves toward the minimum at the origin. This type of model (with very 
small j3) was constructed to account for the possible evolution of the fine structure constant. 
In panel b), we show the corresponding evolution when (3 takes values 0.1, 1, and 10. We 
see that as (3 is increased, the evolution of begins at higher redshift. When (3 = 1(10), 
we also see that the field undergoes a few (many) oscillations about the origin. 




FIG. 3: Evolution of 99* as a function of the redshift for models with potential and coupling defined 
in Eqs. ()3.4ti3.5p . (Left): we set /3=0 and A = 1, 2, 3, 4, 5 and 6, so that we are dealing with a 
quintessence model and no modification of general relativity. (Right): we set /?= 0.1, 1. and 10. 
without a potential (solid) or with a potential and A = 6 (dashed). Differences with and without 
a potential appear only at the very late times. The vertical dashed line corresponds to the time of 
n/p freeze-out and the shaded region to the epoch of BBN. 



IV. NON-UNIVERSAL MASSLESS DILATON WITH QUADRATIC COU- 
PLINGS 

We now consider models in which Ay ^ Ad. For simplicity we focus on models with 
a constant potential and we assume that the two coupling functions are quadratic, which 
gives us a natural extension of the model studied in Ref. [22]. We set 

Av = exp(^^Pv^'?j, AD=exp(^^Pn^'?j, V = Vo, (4.1) 

and we define 

ai{^^) = ln[Ai{^^)]. (4.2) 
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This is also a generalization of Ref. [32[ which was restricted to Brans-Dicke models. 

Vo is constant so that the scalar degree of freedom remains massless. As in the previous 
model considered [Eq. (13. 4j it is associated with a constant energy density in the Einstein 
frame 

PA, = V^/A-kG, , (4.3) 

but scales as Ay'^ in the matter Jordan frame. It follows that the total energy density is 
Pt = P + Pa- Note that the couplings are such that this component is not coupled to ip^, 
that is 

= 1 , «A = 0. (4.4) 

We note that the tests on the deviations from general relativity in the Solar system involve 
only ay and /Sy since they are performed with ordinary matter. 

A. Dynamics 

1. Reduced equation 

The dynamics of this non-universal system can be discussed using p as a time variable. 
Since dp = if*dt*, we can deduce that ip* = H^:ip'^ and we can rewrite the Friedmann 
equation (1A3P as 

(3 - if'^)H^, = S-kG^Xp, + PA*). (4.5) 
The derivative of Eq. (]A2p with respect to p implies that 

^* [if ^^J'^*^ — IP' 1^ ^'^^ " 3w;i)p„. 
Now, it can be checked that 

f = -^'(3-y,:2)(l+*)v.:-^:^ 

where w includes the contribution of the cosmo logical constant, 

p,(l +w) = (p, +PA*)(l + w) . (4.6) 
Thus, we conclude that the dynamics is described by 

-^y.: + (1 - u;)y.: = -5^a.(l-3^.)^^ (4.7) 
3- ip'^ ^ PA* 

= -av-^^^^ - anil - 3wd)-^^ (4.8) 
p* + Pa* P* + Pa* 



This generalizes the equation used in Refs. |16l . |21| . |22| to the case of two couplings. We 



emphasize that this equation holds only when = Vq is constant and for K = 0. 

It is interesting to note that if the cosmological constant is negligible and if the dark 
matter and the baryonic matter have the same equation of state wd = w then Eq. (14. 7p is 
equivalent to the one that would be obtained with one single fluid with equation of state w 
with a universal scalar-tensor theory with an effective coupling function 

A^s{ip,) = Aviip,) + EoADiif.), (4.9) 

which depends on the ratio of dark matter and baryonic matter, as well as initial conditions, 
through Eq. 
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2. Phase space 



Neglecting the cosmological constant, which dominates only in the last e-fold, Eq. (14 .yp 
reduces to 

2 ft,, _l_ /9„'=:„pi(fe-/3v){</'^-¥'^o) 
-(p^ + [1 - w)^^ = — _ ^ , , TT'^*- ^-^"^ 

3 — (f[f 1 + ~^^Q2yf^D-Pv){'Pi-'fit,o) _)_ £iad 

pmat 

Let us focus on the matter dominated era. Eq. (14.101) indicates that the attraction toward 
general relativity will depend crucially on the relative signs of jSy and Pd and on the form of 
the effective coupling function A^fi which enters on the r.h.s. of the Klein-Gordon equation. 
We have the following different possibilities: 

• /3y = 0: the scalar-tensor theory is pure general relativity in the visible sector. Only 
cosmology can set constraints on jSo while all Solar system constraints are always 
satisfied. 

• /?D = 0: the theory is attracted toward general relativity only if /3y > 0. 

• (3d > and /?y > 0: ay and have the same minimum at ip^ = so that the 
effective function aes also has a unique minimum in (p^, = which is an attractor of 
the dynamics. The scalar-tensor theory is thus attracted towards general relativity, 
both for the matter and the dark sectors, during electron-positron annihilation and 
the late time matter era when if^^^O. 

• /?D < and /3y < 0: Neither ay nor ao have a minimum and thus aes does not have 
a minimum, y?* runs to infinity when either matter or dark matter is dominant. The 
scalar-tensor theory thus drifts away from general relativity. We shall thus discard this 
case since local constraints can be satisfied only at the price of an extreme fine-tuning 
of the initial conditions. 

• PdPv < 0: This situation is more complex. = is always an extremum of a^s but 
not necessary a minimum. Two conditions have to be considered. 

First, if 

t-^^ln(--^)<0, (4.11) 

Oeff has a unique extremum at ^9^, = otherwise it has two others given by 

■^....='^o. + ^^^ln(-^— j. (4.12) 

Second, if 

/3 = /3y + So/3d > 0, (4.13) 
the extremum at ^9^, = is a minimum. 

The signs of these two functions set the shape of the function Oeff which determines 
the late time attractor of the Klein-Gordon equation. Note that $ changes sign when 
(5 does. Fig. H] summarizes the various possibilities in the phase space while Fig. [5] 
depicts the modification to Oefr when the values of [3d and /3y vary and when the signs 
of $ and j3 change. 

Given the shape of agg, we have two possibilities for the late time dynamics: 
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Late time attraction to GR j-^ 




FIG. 4: The {Pv,Pd) parameter space with the signs of $ and (3 [defined in Eqs. ()4.1ip and 
(j4.13p ] which determine the shape of agg and thus the late time attractors of the Klein-Gordon 
equation. Three possibilities appear: late time attraction toward general relativity for (3 > 0, late 
time attraction toward a Brans-Dicke theory when f3 < and (3dPv < and runaway when [3d < 
and (3v < 0. The field is respectively attracted toward 0, ±(^att and ±oo. 



j8d = 1, j8„=-1,-5,-10 y8v=10, /3d=-1,-3,-5 




FIG. 5: Modification of the shape of a^s with the values of and fiy and with the signs of ^ and 
p. 

— v?* = is the only minimum. The Klein-Gordon equation involves A^s so that 93* 
is attracted toward during the matter era. Since this corresponds to the mini- 
mum of ay and an, the scalar-tensor is also attracted toward general relativity. 

— (f^: = is a local maximum. In this case (p^ is attracted toward <^att 7^ 0. This 
does not coincide with a minimum of ay and ay — ^ fivfatt- The scalar-tensor is 
thus attracted towards a Brans-Dicke theory, that is a scalar tensor theory with 
lnA = X(f^. Fig. [6] illustrates such a case where y?* cannot relax to the minimum 
of the coupling function Ay because the dynamics is dictated by A^fi. It follows 
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from Eq. fHTTD that 



According to the value of the parameters, the asymptotic Brans-Dicke parameter 
can be compatible with the observational constraints. Indeed, the theory can be 
temporarily attracted toward general relativity due to the mass thresholds in the 
primordial universe if Pv > 0. In this case, ayo can be low enough to pass Solar 
system tests. 

2 . 5 

2 

m 1.5 
01 

tri 1 

Q 0.5 

; ° 
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FIG. 6: When the minimum of Ocff (dashed line) differs from the one of ay (solid line) and 
(long dashed line), the dynamics of the scalar field drives to a minimum of ^eff so that the 
scalar-tensor theory is attracted toward a Brans-Dicke theory. 

As discussed above. Figures H] and [7] summarize the field dynamics. Three possibilities 
appear: late time attraction t_owards general relativity for (3 > 0, late time attraction toward 
a Brans-Dicke theory when jS < and (3dPv < and runaway when (3d < and jSy < 0. 
In each of these cases, the field is attracte_d towards 0, iy^att and ±oo, respectively. Only 
models such that the effective parameter (3 is positive, enjoy a late time attraction toward 
general relativity. 

Fig. [7] connects the phase space {I3v-,(3d) with specific examples of model evolution which 
highlight the different possibilities. Fig. H] has been miniaturized in the lower left of Fig. [7] 
and with evolutionary examples for each quadrant except that for which Pd, f3v < and the 
theory exhibits runaway behavior. 



In the upper left of Fig. [TJ we show the evolution of ipl as a function of redshift for 
P = —1, 0, +1, all with Pv = —10. If Pv is negative the theory is dragged away from general 
relativity. This is seen by the increase in ipl prior to BBN. This could be compensated at 
late times if /3 > as seen in the lower curve with P = 1 which undergoes two oscillations 
about the minimum. Note that when P = the Klein-Gordon equation has no source term 
even after the end of the radiation era so that y?* remains frozen to the same value than 
during the radiation era. The scalar-tensor theory can be attracted toward general relativity 
only due to the effects of the mass thresholds. Therefore no evolution is seen after BBN 
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FIG. 7: (Lower left): Phase space of the solutions. (Upper-left) when /3y < the theory is dragged 
away from general relativity during BBN and then, depending on the sign of /3, attracted towards 
or away from general relativity at late times. The field is at the minimum of the coupling function 
today only if /? = (as seen by the horizontal line at late times). (Lower-right) when (5d < 
the theory is attracted toward general relativity during BBN and then, according to the sign of 
f3 attracted or not toward general relativity at late time. (Upper-right) Both /3s are positive and 
the theory is always attracted toward general relativity. The opposite is true when both (3s are 
negative. 



when j3 = 0. When (3 = —1, the field continues to evolve away from ipl = when matter 
domination sets in. However, and as discussed in detail in Ref. |2^, there are many other 
mass thresholds prior to electron-positron annihilation and a. jSy < model would have to 
be fine-tuned to be attracted close enough to general relativity today. We shall thus restrict 
our analysis to /^y > 0. 

Models with f3v > are attracted toward general relativity during electron-positron 
annihilation (see discussion below). Then, if (3^ < 0, the theory is either attracted toward 
general relativity or a Brans-Dicke theory depending on the sign of (3. As the analysis above 
also shows, in the latter case the theory is attracted toward a value of (p^ different from 
(p^ = 0. Examples of models with Pv > 0, but Pd < and /3 = —3, 0, -|-3 are shown in the 
lower right plot of Fig. [71 Each case moves towards the origin (and undergoes one oscillation 
about the origin) prior to and during BBN. At late times, we see again that the model with 
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/3 > continues to evolve towards the origin after matter domination, while model with 
(3 < evolves away from general relativity. As before the case with j3 = shows no further 
evolution. 

Finally, in the upper right panel of Fig. [3, we show an example of a model with (3d = 
j3v = +10. The field moves through the origin once during electron-positron annihilation, 
and then later continues to oscillate about the origin. For completeness, we also show an 
example with = Pv = —10 in this panel which shows the run away behavior at late 
times. 

In conclusion, we will concentrate our analysis on the case in which both jSy and (3d 
are positive (see upper-right plot of Fig. [T]). This corresponds to models that are attracted 
toward general relativity without involving a tuning of the different parameters and the 
initial conditions for the value of (3 to be positive. 



3. Different regimes 

Deep in the radiation era, the cosmological constant the dark matter and baryonic matter 
components are negligible so that Eq. (14.71) simplifies to 

^ + = 0. (4.14) 



As a consequence, the field is frozen at a constant value and the initial conditions can be 
chosen as 

^'*m = 5 'f* = = constant . 

Still in the radiation era, when the universe cools below the mass of some species, Xi this 
species becomes non-relativistic and induces a non-vanishing contribution to the r.h.s. of 
Eq. (14.71) . The most important mass threshold for the BBN predictions is the last of these 
thresholds, the one associated with electron-positron annihilation. Previous mass thresholds 
are also important and drive the scalar-tensor theory toward general relativity (see § III. A. 3 



of Ref. |22|). We assume that the dark sector is not affected by any mass thresholds after 
muon annihilations. In that case, we can neglect the effect of the dark sector and we end 
up with the same result as in Ref. j22|, 

^ + ?<^: + S,(T)/5y(^, = , (4.15) 



■ 3 

with 



This implies that we are driven toward general relativity during the electron-positron anni- 
hilation only if /5y > (see right half of Fig. [Tj). 

After electron-positron annihilation, the universe is still dominated by radiation and the 
scalar field freezes at a constant value, v'*out- BBN can place a constraint on the value 
of ttout = o(v^*out)- Unfortunately, this constraint depends on Oin which is unknown. To 
compare these constraints to those obtained in the Solar system, we need to relate Oout to 
ao- Thus, our code integrates the evolution equation up to the present, so that we obtain 
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Oo directly. For the particular case of a vanishing potential or as long as the field is slow 
rolling, y)' <^ 3, Eq. (14. 7p takes the slightly simplified form 



2 /3y + /3z)Soe^(fe-^^)('^>-^*o) 

o<^* + (1 - w)(p^ = — - — _ 1 — ^1?:;;^^:;^^'^* 



(4.17) 



Contrary to the case of a universal coupling (§ III. 5 of Ref. [2j]) this equation cannot be 
integrated analytically in the matter-radiation era. 

The field is attracted toward its minimum and at late times we can assume that ^ 1, 
in which case the equation of evolution reduces to 



yi + (1 - w)ip', + ^v* 



0, 



2 (^V + Pd^O 
UJ = = 



It follows that the field evolves as 



e 4> 



e 4' 



, 3 \ „ /3 
A cos ( -rp j + I -rp 



A cosh ( ^f'Pj + -Bsinh (^^P 



respectively if cj^ < | and cj^ > | with 



This allows us to compute the period of the last oscillations, 



Ap= - 



71 



3r(/?y,/?^,So) 



(4.18) 
(4.19) 



(4.20) 



Unfortunately we can not compute the phase, which is required in order to match to the 
solution in the radiation era. 

An example of such an evolution is depicted in the left panel of Fig. [Hlwith jSy = Pd = +^0 
as in the upper right of Fig. [71 Here we see the same oscillations scaled by e^^^^ as expected 
from Eq. (14.191) with period given by Eq. (14.201) . In the right panel of Fig. [8], we show several 
cases (as labelled) showing the late time scaling and oscillations. We can see that the field 
remains frozen during the radiation era up to the kick during electron-positron annihilation. 
When the universe starts to be matter dominated the field undergoes damped oscillations 
with a period given by Eq. (14.201) . 

B. Constraints 

1. BBN 



We proceed as in Ref. [22| and vary the baryon to photon ratio, r;, the values of (3v and 
a.s well as the initial conditions but we assume that the other cosmological parameters 
are fixed to their standard values. We compute the light element abundances as a function 
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FIG. 8: (Left): Evolution of exp(|p) as a function of redshift. The frequency of the oscillations 
depends on /Jy, (3d and the energy densities as in Eq. (j4.20p . (Right): Some solutions with (3 > 
summarizing the effect of Pv during BBN and the evolution during the matter era. 




FIG. 9: (Left): Light element abundance as a function of ain when rj = t/wmap and Pv = 10. 
(Middle): Same as a function of (3v when r] = r/wMAP and Oin = 0.1. (Right): Same as a function 
of r] for Py = 10 and Oin = 0.1. The standard model is shown by the lower (blue) dashed curves. 

of these four parameters. Fig. [9] shows the effect of the independent variation of each of 
these parameters. This allows one to set the BBN constraints on these parameters which 
are then propagated until today to infer constraints on ayo, Pv and Pd- Fig. [TO] summarizes 
these constraints on the {i3v,Pd) plane obtained from BBN. 

In the left panel of Fig. [9|, we show the resulting light element abundances as a function 
of the parameter for fixed values of Pv = Pd = 10 and t] fixed at the WMAP value of 
6.23 X 10~^° [4]. As one can see in this panel and others to the right, there is very little 
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dependence of ^He, D, or ^Li on either Oin or Pv Note that the ^Li abundance is always 
in excess of the observations for this value of 77 (see e.g. Ref. 48|). In contrast, there is a 
relatively strong dependence of the ^He abundance with 0^. The shaded regions correspond 
to the range of the observational determinations as described in Appendix C. The upper 
limit on the ^He abundance of 0.259 places a constraint on a^^ < 0.2, when f3v = 10. In 
the middle panel of Fig. [HI we show the element abundances as a function of /3y for fixed 
= 0.1 and the WMAP value for 77. At this value of Oin, we find a lower bound on jSy > 2.8. 
Finally, in the right panel of Fig. [9|, we show the abundances as a function of 1] for fixed 
Oin = 0.1 and jSy = (3d = 10. This choice of parameters leads to an increase in the ^He 
abundance of roughly 0.005 over the standard model with Oin = and /5y = as shown by 
the lower set of blue dashed curves. 




FIG. 10: Constraints on ayo in the plane (Pv^Pd), respectively for logayo = —5 (solid line), —4.5 
(long dashed lines) —4 (dashed lines) and —3.5 (dotted lines). The red line in the lower left corner 
shows the solar system constraint of logayo = —2.5 



Now, for each value of the pair {Pv, Pd), we can compute the abundances of "^He, ^He, D 
and ''Li as a function of which allows us to constrain the values of this parameter from the 
observational data. We have set rj = ?7wmap, though in principle, one could derive combined 
constraints on the pair of inputs (ain,?]). One can check that this is a safe approximation 
and that the deviations fro m g eneral relativity for the models considered do not lead to 
significant CMB deviations |23j. Then, this constraint on can be propagated to get a 
constraint on the value of (f^, today, or equivalently on ayo- 

The result of this analysis is depicted in Figure [lOl This figure can be understood 
as follows. First we see that almost everywhere on the plane the solar system bound, 
< 10~^, is satisfied (the area above the red line barely visible in the lower left corner); 
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see Eq. ( ]B2p . Second, for a given value of {i3v,i3d) we can read off the maximum value 
of deviation from general relativity that can be achieved in the Solar system once BBN 
constraints are satisfied. The result from Ref. [l^l corresponds to the cut (3v = Pd of 
Fig. [TUl Figure [TD] shows clear structure with two sets of minima: (1) a set of parallel 
lines with increasing periodicity and (2) a vertical line close to /5y ~ 5. Let us now try to 
understand this behavior. 

First, the periodicity seen in this figure is not directly related to the period fl4.20p of 
the time evolution of ip^, but rather on the phase of this solution when evaluated at p = 0. 
Since we do not have an analytical solution in radiation-matter era, as in Ref. |22j, our 
discussion can only be approximate but will still shed some light on our result. During the 
radiation era, but after BBN, the field is frozen at a value v^out*- Assuming that at the time 
of matter-radiation equality, p^q, (peq* — V'out*, we can obtain a solution for ip^ at late times 
by matching the solution (I4.18p to the constant solution at Peq so that 



-f (P-Poq) 



COS 



4 



It follows, then, that the value of the scalar field today is given by 

'3 



^0* 



^ouuG^^"" COS 



rrpeq 



(4.21) 



(4.22) 



which illustrates that the oscillations of (po^, and thus of avo, as a function of the two 
parameters {Pv,Pd), are due to the oscillations of ip^^o as a function of r and are thus 
determined by the phase at the time of matching. To determine the points at which ip^,o = 0, 
we need the value of peq that is determined by imposing that pmat + Pd = Prad- This leads 
to 

The almost horizontal periodic structure in Fig. [10] is related to the solution such that 
V^o* = 0. When this happens, Eq = Eq and ^^(v^o*) = and these zeros occur when 
nvr with n eN, that is when 



8 f3v + PdEq 



3 1 + 



- 1 



TT 



-n- 



+ oOvout + In 1 + Soe i^v 



fip-fiy 



Now, since the dark sector does not influence the dynamics during electron-positron annihi- 
lation, (y9out* depends only on the values of ipin* and This complicated relation is depicted 
in Fig. 2 of Ref. [23]. From this analysis, we argued that ayout is at most of order 10~^-10~^. 
This implies that we expect the denominator of the r.h.s. of the previous expression to be 
dominated by In . Thus, the positions in the plane Pv^ Pd for which ayo = (and as a 
consequence aco = 0) corresponds to a set of parallel lines C„ determined by 



Pv + PdEo = ^(1 + Ho) 



Ann 



31n- 



+ 1 



(4.24) 



This explains why the spacing between C„ and Cn+i grows almost linearly and is mainly 
determined by Eq and f^mato/^radO) which arises from the fact that the oscillations are related 
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by the time of matching. The slope is determined by Sq since it characterizes the evolution 
during the matter era. 

Let us now turn to the second feature in Fig. [TOl that is, the vertical line on which 
avo = 0. The relation between ipo^t and v?in is not monotonous and is determined mainly 
by j3v (see Figure 2 of Ref. 22(] that shows that it is periodic in (3y with a first minimum 
for Py ~ 5). Figure [TT] shows how that the relation between 99 out and ipiri has indeed a 
minimum at /3y ~ 5. This implies that for these models lp^, freezes to a very small value 
after the electron-positron annihilation so that these models are more efficiently attracted 
toward general relativity. 

In summary, we can obtain a good understanding of our numerical results and of the 
origin of the structures of the minima in Fig. [TOl In particular, we found that the two 
sets of structures do not have the same physical origin, one is related to electron-positron 
annihilation and the second to oscillations in the matter dominated era. 



0.03, 0.1, 0.3 




FIG. 11: Evolution of (font* as a function of /3y assuming that = 0.03, 0.1 and 0.3. We see that 
the first minimum is obtained for f3v ~ 5, which corresponds to the vertical minimum of Fig. [TOl 
There is a slight dependence of the minimum with respect to (pi^ and hence with the constraints 
provided by BBN. 



2. Time variation of the gravitational constant 

From Eq. (]B3p . using Eq. (I2.19p . the time variation of the gravitational constant is given 

by 

a = 2ay (1 + (4,25) 

where a = if'^dlnGcav/dt is constrained by Eq. (IBSp . i.e. we must have |cro| < Sq, where 
So is the experimental upper bound on a today; see Eq. (IBSp . 
Assuming we are in slow-roll today, we deduce that 

(1 + l^Ao)<^*o VO- 

-L + An 
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We also deduce that the time variation of the gravitational constant today is 



1 + 



1 + Ho)(l + ^Ao) - a^^iPv + Pd^o 



(4.26) 



Figure [T2] shows the region of the parameter space for which |(Jo| < So given the 

deviation from general relativity today, i.e. ayo- The contour plots can be understood as 
follows. Since Oyg < 10~^, the second term in the denominator of Eq. (14.261) is negligible. 
Unless we have /3s of order 10^, we conclude that 



1.6 X 10-"ayo(l + Pv){(3v + Pd^o)- 



This shows that even if the dynamics of the field is influenced hj jSn which can be very large, 
the Solar system constraints on ayo are so strong that it requires very high to obtain a 
sizable time variation of the gravitational constant. 




FIG. 12: Contour plot of the time variation of the gravitational constant in the plane (fiv^f^D) for 
different values of log ayo indicated on each curve. The model has to lie on the lower-left part of 
the plot so that \ao\ < Sq. 



V. CONCLUSIONS 

For a theory of gravity universally coupled to matter, a scalar-tensor theory can be defined 
by a single coupling function A{(f) and a scalar potential. Such a theory can affect the 
evolutionary history of the universe, particularly in models of extended quintessence. Here 
we examined two types of quintessence models and examined the effects of the non-minimal 
coupling on the evolution of the quintessence field. For inverse power law quintessence 
models, we showed that the evolution towards a tracking solution is accelerated due to mass 
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thresholds in the radiation dominated epoch, and during the subsequent matter dominated 
epoch. In contrast, we saw that for quintessence models with defined minimum which is 
also a minimum of the coupling function, the non-minimal coupling to gravity has very little 
effect except at late times. 

It is, however, possible that the dark sector of the theory couples to gravity differently 
from that of the visible sector. In this case, there is an additional coupling function as 
defined in Eq. fl2.1l) . In this article, we have worked out the cosmological evolution of a 
scalar-tensor theory with distinct coupling functions for the coupling of the visible and dark 
sectors of the theory to gravity. We developed the qualitative features of the evolution 
with respect to the parameter space defined by /3y and Pd which are derived directly from 
the two coupling functions. We delineated those regions of the parameter space which are 
attracted towards general relativity, towards Brans-Dicke gravity, or towards a runaway 
solution. For those theories with a late time attraction to general relativity, we also derived 
additional constraints on the parameter space from BBN and precision gravitational tests. 
In particular, it was shown that in this class of models constraints more stringent than those 
derived from the Solar system can be obtained from BBN. This provides some insight on 
the equivalence principle between the visible and dark sectors. 
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APPENDIX A: NUMERICAL IMPLEMENTATION 

From a numerical point of view, it is easier to integrate the Einstein equations in the 
Einstein frame while using the cosmic time t in the MJF. The expressions for the energy 
densities of the visible sector are trivial in the MJF, where they can be integrated analytically. 
It follows that the system reduces to 
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_ 

dt 

Hl = 

dprad 

dt 
dpmat 
dt 

+ 3ifp£, = Ay^ {an - ay) PdV'*, 



^ 1 2 2 

+ Prad) + + 3^ ~ 



matPmat J 

K 



+ 4iyprad = 
+ 3iJpmat 



(Al) 
(A2) 

(A3) 

(A4) 

(A5) 

(A6) 
(A7) 



Note the asymmetry in the couphng functions, A, arises from the fact that we have inserted 
the MJF energy densities rather than those defined in the EF. The last equation is necessary 
to derive the Hubble parameter in the MJF where the nuclear reactions are integrated. 



APPENDIX B: LOCAL CONSTRAINTS 

The deviation from general relativity are constrained in the Solar system. These devia- 
tions are usually summarized by constraints on the post-Newtonian (PPN) formalism [49|]. It 
is a general formalism that introduces 10 phenomenological parameters to describe any pos- 
sible deviation from general relativity at the first post-Newtonian order j49|. The formalism 
assumes that gravity is described by a metric and that it does not involve any characteristic 
scale. In our particular case, it is necessary that the scalar field is light so that a Yukawa 
interaction on the scale of the solar system is not induced. 

Since all of the matter in our Solar system stems from the visible sector (i.e. the dark 
matter component of the Sun and planets is supposed to be negligible), the two Eddington 



parameters (see Refs. [49|, l50|) can be expressed in terms of the values of ay and (3v today 



as 

Solar System experiments set strong limits on these parameters. The perihelion shift of 
Mercury implies j5l| \2'-^^™ — j3^^^ — 1 1 < 3 x 10~^, the Lunar Laser Ranging experiment 13 
sets 47^^^ — 0^^^ — 3 = —(0.7 ± 1) X 10"'^. Two experiments give a bound on 7^^^^ alone, 
the Very Long Baseline Interferometer j53| sets 17^^"^ — 1| < 4 x 10~^ and the measurement 
of the time delay variation to the Cassini spacecraft near Solar conjunction 5J] 7^^^ — 1 = 
(2.1 ±2.3) X 10"^ 

These two last bounds imply that ayo has to be very small, typically 

ayo < 10-' (B2) 



while /Syo can still be large [55j. Binary pulsar observations imply /5yo ^ —4.5. Note that 



even though (3o is not bounded above by experiment, we will assume that it is not very large, 
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typically we assume /3vo ^ 100, so that the post-Newtonian approximation scheme makes 
sense. Note that none of these observations constrain the dark sector and (azpoj/^Do) are 
completely free. 

Since Cavendish experiment and planets are composed only of baryonic matter, we deduce 
that the gravitational constant is given by 



Gr 



G,A'y{l + a'y] 



(B3) 



Using p as the time variable as defined in Eq. (12.201) . the time variation of the gravitational 
constant today determines the time derivative of y?* today 
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dp 



Experimentally, this is bounded 5^ by 

1 dGnav 



Gr 



dt 



iSnl < 5.86 X 10-^h-\ 



(B4) 



(B5) 



The new main feature compared to the models analyzed in Ref. [22| is that the dynamics of 
the scalar field depends on the dark sector coupling and thus on a^- 



APPENDIX C: BBN CONSTRAINTS 



BBN is one of the most sensitive available probes of the very early Universe and of physics 
beyond the standard model. Its success rests on the concordance between the observational 
determinations of the light element abundances of D, '^He, ^He, and ^Li, and their theoreti- 



cally predicted abundances |57l.l58l|. Furthermore, measurements of the CMB anisotropies by 
WMAP j^, 0] have led to precision determinations of the baryon density or equivalently the 
baryon-to-photon ratio, rj. The new WMAP 5-year data alone is Q]^h'^ = 0.02273±0. 00062 
and is equivalent to r/io.cMB = 6.23±0.17, where rhn = 10^°?7. Using the WMAP data to fix 
the baryon density, the light element abundances H, 6^, 61, 62, 63] can be quite accurately 
predicted. The wealth of the cosmological data are obtained from spectroscopic observations 
and compared directly with BBN predictions assuming the WMAP determination of Q^^h"^. 

Note that the "^He abundance is often used as a sensitive probe of new physics ( see e.g. 
(63|). This is due to the fact that nearly all available neutrons at the time of BBN end up in 
^He and the neutron-to-proton ratio is very sensitive to the competition between the weak 
interaction rate and the expansion rate. 

In the following, we briefly state the abundance measurements used in our analysis. 



a. D/H 

The best determinations of primordial D/H are based on high-resolution spectra in high- 
redshift, low-metallicity quasar absorption systems (QAS), via its isotope-shifted Lyman-o; 
absorption. The seven most precise observations of deuterium in QAS give (see [65| and 
references therein) a weighted mean value of D/H = (2.82 ± 0.21) x 10~^ (l^cr) good 
agreement with BBN at the WMAP value of t]. Note that the uncertainty quoted above is 
purely statistical and there remains considerable scatter in the data. 
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b. ^He 



^He is observed in clouds of ionized hydrogen (HII regions), the most metal- poo r of which 
are in dwarf galaxies. There is a large body of data on ^He in these systems 66,l67l| for which 

However, 



m 

an extended data set including 89 HII regions obtained ¥„ = 0.2429 ± 0.0009 
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the recommended value is based on the much smaller subset of 7 HII regions, finding Yp = 
0.2421 ± 0.0021. 

^He abundance determinations depend on several physical parameters associated with 
the HII region in addition to the overall intensity of the He emission line. These include, 
the temperature, electron density, optical depth and degree of underlying absorption. Un- 
fortunately, there are severe degeneracies inherent in the in the determination of the ^He 
abundance [6^. Using a subset of the highest quahty data from the sample of Izotov and 



Thuan [66|, Monte Carlo methods were used to extrapolate the ^He abundance 69|] which 



was determined to be 1^ = 0.2495 ±0.0092. Conservatively, it would be difficult at this time 
to exclude any value of ¥„ inside the range 0.232 - 0.258. 



c. ^Li/H 

The systems best suited for Li observations are metal-poor halo stars in our Galaxy. 
Analyses of the abundances in these stars yields [70| Li/H|p = (1.23lgjg) x 10^^° when 
systematic uncertainties are included. The ''Li abundance based on the WMAP baryon 
density is predicted to be ^Li/H = 4.15l|^:^^ X 10 We note that a recent reanalysis 
of the '^He(a, 7) ''Be reaction, which is the most important ''Li production process in BBN, 
was considered in detail in [711]. When the new rate is used a higher ^Li abundance is found 
H ^Li/H = (5.24^°:^^) X 10-10. 
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